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Abstract 

We study dynamical properties of bouncing particles inside of channels with 
sinusoidal walls. Taking as parameters the amplitude and the phases between the 
walls we study the transmitivity and its dependence on these parameters. We 
find an analytical approximation for transmitivity for small amplitudes, which is 
corroborated by numerical calculation. 
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1 Introduction 



Two dimensional billiards are systems consisting of a free point particle moving in a 
two dimensional region bounded by rigid walls. The billiards problem leads to both 
classical and quantum interesting and well studied class of Hamiltonian systems such as 
the stadium [1] and the Sinai billiards [2 j. Basically, from a classical point of view the 
problem consists in to find the trajectory of a particle which bounces elastically on the 
boundary walls. From the quantum point of view, the problem consists in solving the 
Helmholtz equation by requiring that the wave function vanishes on the boundary. At the 
mesoscopic level, a question of great interest in the transport properties in open billiards, 
such as quantum dots is the influence of the shape of the boundary. In fact, it has been 
observed that coherent backscattering and conductance fluctuations are strongly shape 
dependent jl] . Due to the development of nanostructure devices El (in addition to the 
quantum dots we can mention the quantum wires), the transport properties with both 
regular and irregular profiles is of current interest |H1 El • As matter of fact it has been 
observed that when the nanotubes are elongated they present a pattern of regularity in the 
shape. Many of these systems have been studied and show regular and chaotic dynamics 
jH]. Quantum and classical billiards systems, have been helpful in the understanding 
the quantum manifestation of classical chaos and quantum transport for large values of 
energy. At low temperatures, of order of few mK, the mean free path of electrons is much 
larger than the size of the quantum devices, diffusive effects disappear and the electrons 
show a rather ballistic behavior. In this case a classical analysis is enough for studying 
the transport properties in these quantum and mesoscopic devices. 

In this work we study a two dimensional rippled open channel. The profile of walls of 
the channel are sine-like but, in general, with different phases. The case of one rippled wall 
has been previously studied in Ref. [EllHj- We investigate the dependence of the transport 
properties on the difference of phases as well as the dependence on the amplitude between 
the walls in the ballistic regime, we obtain an analytical approximation for reflexivity for 
small ripples. We observe that the chaotic regime appears faster when the difference of 
the phase between the walls is near to tt than when they are near to be in phase for a 
given amplitude. Our results are corroborated by numerical calculations. 

2 The rippled channel 

Let us consider an open channel with rough walls such as the shown in Fig. ^ The upper 
and lower walls are determined, respectively, by 

yi = b + a sin 2ttx 

2/2 = — b + asin27r (x + r) (1) 

where, instead of the physical size, we use dimensionless variables x = y,y=j,6=j, 
a = j, L = H- with / being the length of the period and L' the length of the channel, 
A, B are respectively the separation between the walls and the amplitude. The variable r 



2 



denotes the difference of phase between the upper and lower walls and it is restricted to 
the values | r |< ~. 

The trajectory of each particle bouncing elastically inside the channel is define by the 
initial condition (y , ao), where yo denotes the high at which the particle was initially 
dropped and «o the initial angle. Since the particles bounces elastically, it satisfies that 
the angle of incidence equals to the reflection angle, meseared from the tangent line to 
the curve. For our work an important transport property is the transmit ivity, T, that is 
defined by the flux of the transmitted particles divided by the incoming flux. To compute 
this numerically we dropped 10 4 particles from the left side of the channel (at x = |), of 
length L = 2, by 10 sources uniformly distributed along the y axis (each source eject 1000 
particles). We assume that the particles do not interact among themselves and that no 
multiple collision with the walls occurs before the particle leaves away the corresponding 
wall. 

The angular distribution (for each source) of the dropped particles is 

p("o) = ||cos(a ), (2) 

where ao is measured with respect to x axis as is shown in Fig. ^ and uq is the total 
number of particles. 



3 Poincare sections 

Let us discuss a channel whose walls are separated away by a distance 2b = 0.2. In 
order to obtain all the different types of trajectories in the phase space (x n ,p n ), were 
x n and p n are respectively the value of the x coordinate at the n-th collision with the 
upper wall and p n the x component of the momentum after the n-th collision. We throw 
particles with different initial conditions in the phase space. The x component of the 
momentum is normalized by the velocity. The system can be carry out in practice if we 
identify the particles as electrons and resort to the ballistic regime which is realized at low 
temperatures in nanostructures where the size of the sample L' is such that A e << L' < L e , 
where A e is the electron wavelength and L e is the mean free path 6J. 

Some Poincare sections are shown in Figs. 0E Notice that Poincare sections are 
moved along the x direction by a quantity j. For small values of a the corresponding 
Poincare sections are similar to a pendulum phase space (Figs. Eland EI). When we go to 
higher amplitudes, the separatrix in x n ~ 0.5 of Fig. |3]becomes chaotic. For higher values 
of a we get a region of trapped particles (ellipses or deformed ellipses) surrounded by a 
chaotic sea (Fig. |7J), this occurs when the last KAM curve breaks and then all chaotic 
regions are connected. For small values of a, the elliptic orbit indicates that the particles 
are trapped around the center of the ellipse (not all of these points are fixed points). 
Notice that for a fixed point Y n = Y n+ i, where Y n is the trajectory of the particle after 
the n-th collition. This leads to the condition X = ^Arcos (4^); where X is a fixed 
point, so in order to have a fixed point the phase r must satisfies the relation |r| < Airba. 
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When the amplitude is small, the trapped particles execute librational motion (similar 
to the pendulum oscillations around the stable equilibrium point), they move adiabatically 
forward and backward. The only change due to the phase, for small amplitudes, is that 
the center of the ellipses are moved in the x direction but they maintain almost constant 
in the p direction, near to p n = (see Figs. EJandlHJ). 

From our results, we observe that for values of the phase near to r = |, whenever 
the amplitude is increased, the chaotic regime appears faster than for values near to 
r = 0. This behavior is observed as well for the last KAM curves; as we go near to 
r = | the last KAM curves break sooner. As an example when the walls are in phase the 
separatrix becomes chaotic until a = 0.04 and the last KAM curve breaks until a = 0.075. 
For r = 1/3 the last KAM curve breaks around a = 0.016 and for greater values, for 
instance a = 0.02 (FigEJ) it appears a cluster of seven elliptical curves surrounding the 
principal first order resonance island. When r = 1/2 the chaotic behavior appears around 
a = 0.005 and the last KAM curve breaks around a = 0.015. There is an appreciable 
difference between region of trapped particles when the profiles are in phase than when 
r = 0.5. For the same amplitude a, we can see these differences from Figs. El and 0] there 
are more trapped particles when the phase is r = 0.5 than when the phase is r = 0. 



4 Transport properties for small amplitudes 

For small amplitudes (e.g. a = 0.001) and when the direction, at which the particles are 
dropped, is almost parallel to the vertical, the particles collide almost perpendicularly 
with the walls, this means that \y\ » \x\. As the particles bounces elastically the energy 
is conserved. If, in addition, we use the principle of adiabatic invariance we can show that 
the motion in the x-direction is described approximately by [10J: 

( x f = — - (-^— V , (3) 
m \D(x) J ' 

where C is a constant. In this case C = yD(x) is the action, E is the total energy and 
D(x) is the separation betwen the walls at the x point, v « y. The term V(x) = (j^y) 
can be interpreted as an effective potential. If | << 1, V(x) becomes 

V(x) w ( X _ | (sen(2vrx) - sen(2vr(x + r))) \ . (4) 

For the case of r = |, V(x) takes the form of pendulum-like potential and explain why 
the phase space in this instance has a pendulum-like form (see Fig. |BJ). We may label, 
in Eq. the variable x with the index n to indicate the ra-th collition at the x point. 
Let x n = vcos(a n ) and x m = vcos(a m ) be the x velocity at the n-th collition and m-th 
collition respectively. From Eq. © we can obtain a relation between a n and a m : 

, 2 (cosV - cos 2 « m ) = (^jL) - , (5) 
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where D n (x) = D(x n ). Let be (3 n = | — a n , the angle that the trajectory makes with 
the vertical at the n-th collition so that (3 n is small and it becomes small and small as it 
reaches the turning point x^ (see Fig. |2J), at this point (3^ ~ 0. Let be (3 C the critical 
angle at which the particle was dropped in x = | and does not turn back. The trajectory 
corresponding to this initial condition is associated with the the largest amplitude of 
librational motion at the center of the ellipses. Let us define x^ m as the point (being 
bigger than x = |) for which the greatest ellipse crosses the x axis. 

Let (3 n = p c , D n = D c = D(\) and f3 m = @n, D m = D(xNm) substituting these into 
Eq. (j3J) and keeping only terms of the order of a and (3% we obtain: 



{(3 C ) 2 ~ 7(sen(27r(xAr m + r)) - sen(27rajv m ) + sen(27rr)). (6) 
o 



If r = ~, then XNm ~ 0.75 and: 

ft 2 « f • (7) 

This result is similar to the semiplane case, the only difference is that in this case the 
amplitude of the channel is twice bigger than the semiplane case jHj. Notice also for r = 
and x^m = 0.75 and from Eq. (jHJ) that (3 C is of order of zero. Here, we can see why the 
great difference between r = and r = 0.5, even for small amplitudes. If r = 4, it can be 
seen from Fig. EJthat x^ m = | and then (3 2 c + l)f , this means that for the case of 

r = \ the channel has a higher value of resistivity than for case of r = |. 

For small amplitudes, we see numerically that (3 C is almost the same for different values 
of yo where yo, is the initial point for which the particle was dropped, this means that 
(3 C depends mainly on r, b and a. We have seen numerically that for small amplitudes 
the phase space is dominated by KAM tori in a pendulum-like manner, the motion of the 
particles with angles ao ~ \ can be considered in the adiabatic approximation, thus the 
velocity in the x direction v x , can be estimated as v x oc (3q and produce a backward flux, 
such that the reflexivity, R, is given 
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Using the approximation for f3 c given in (jSJ), Eq. (jHJ) becomes 

R K h& 



G(r)as, (9) 



where the coefficient 
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G{r) = 3- (sen(27r(a;Ar m + r)) — sen(27TXAr m ) + sen(27rr)) 2 . (10) 
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Notice also from equation Q that for all phases the reflexivity is proportional to a 2 
and the information of the phase is contained in the coefficient G because we have seen 
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numerically that for small amplitudes x^m varies very slow when a is varied, this means 
that x Nrn depends mainly on r. An important measurable quantity is the resistivity p R 
of the channel that is related with the reflexivity through Landauer's formula JT] . 

R 

Poc- (11) 

Since T — > 1 as R — > and according to Eqs (jUJ) and (JUJ, we conclude that pr increases 
proportionally to a a independently of the phase r. 



5 Numerical results for reflexivity and transmitivity 

Fig. |H1 shows reflexivity against amplitude for three different phases r = |, r = \ and 
r = | which correspond to the following numerical equations respectively R\ = 47.46a 1,46 , 

i?2 = 24.53a 1 ' 38 , i?3 = 46.9a 152 . Our analysis predicts R oc at (Eq. |DJ) which gives a good 
agreement between the analytical and numerical results in powers of a. The analytical 
coefficients G(r) for r = \ and r = 7 are respectively 37.9 and 49.7, the discrepancies in 
the coefficients, between the analytical and numerical results arise from the exponential 
grow of the numerical error since we made a fitting data as a function of a power law. 

In Fig. El we show numerical values for resistivity against the phase keeping constant 
the amplitude. We can see that for r = the reflexivity vanishes; a rather expected result 
since from Eq. (JHJ) we obtain that (3 C ~ 0. Moreover when the phase (r) is increased the 
reflexivity increases until it reaches its maximun value at approximately r = 0.35, then 
the reflexivity goes down to the value at r = 1, which is the same geometrical situation 
that for r = 0. Notice that there were significant change in reflexivity only when a change 
on the phase r is more than 8r = 0.05. 

We have seen numerically that T = 1 when r = if a < 0.005, this does not mean that 
all the particles are transmitted, but the critical angle is so small that makes the computer 
insensitive to this values of angles, this implies (3 C < 5 x 10~ 3 , using Eq. Q we see that 
R < 1.25 x 10~ 7 this is in agreement with Eq (jUJ) that predicts B c w 0, which implies 
R ~ 0. We can from Eq. (fTTj) that pn < 10~ 7 . In general the numerical calculations are 
in agreement with the analytical results. 

6 Conclusions 

We have analyzed a two dimensional channel with rough walls. We studied the Poincare 
sections for different values of phase between the walls r and the amplitude a. We found 
that Poincare section are pendulum-like for small values of a. We also note that when 
the amplitude is increased the chaotic regime appears sooner when we are near to r = 5 
than when we are near to r = 0. Based on analytical and numerical results we found that 

3 

resistivity p oc aa is independent of the parameter r, such that the information about 
the phase is complete contained in the proportionality coefficient G(r). Our numerical 
results for small amplitudes, show that the system with the biggest reflexivity is around 
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r = 3.5 and the system with the smaller reflexivity is when the walls are in phase and 
resistivity is of order of a 2 . We have shown that the difference of phases between the 
walls play an important role in dynamic transport properties even for small values of 
the amplitude. Finally we would like to mention that our results could be useful in the 
analysis of transport properties in elongated nanotubes in the ballistic regime. 
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Figure 3: Poincare plot for b = 0.1, a = 0.01 and r = \. 
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Figure 4: Poincare plot for b = 0.1, a = 0.01 and r = 0. 
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Figure 5: Poincare plot for b = 0.1, a = 0.001 and r = \. 
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Figure 6: Poincare plot for b = 0.1, a = 0.001 and r = ^. 
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Figure 9: Numerical values for Reflexivity R against phase r and amplitude a = 0.001. 
Notice that for the values r = and for r = 1 it is same physical situation. 
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